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SECTION  I 
INTRODUCTION 

The  effects  of  external  store  carriage  on  aircraft  performance  can  be 
significant,  especially  for  fighter  aircraft  carrying  fuel  tanks  or  air  to 
ground  weapons.  Degradation  in  top  speed  and  range  must  be  accurately 
estimated  to  determine  if  a  particular  store  configuration  is  practical  from 
a  mission  standpoint.  Flight  testing  can  provide  answers,  but  it  is  extremely 
expensive  and  time  consuming.  For  this  reason,  only  selected  aircraft  loadings 
are  usually  tested. 

To  estimate  performance  parameters  for  a  greater  number  of  configurations, 
wind  tunnel  test  data  are  often  used.  A  scale  model  of  the  clean  aircraft  and  one 
with  the  aircraft  loaded  with  the  store  configuration  of  interest  is  tested. 

Curves  of  lift  versus  drag  coefficient  are  constructed,  and  the  drag  increment 
due  to  external  stores  is  calculated  for  a  particular  flight  condition,  i.e., 
lift  coefficient. 

Because  of  the  nature  of  any  measurement  process,  an  uncertainty  is 
associated  with  the  wind  tunnel  data.  The  total  data  uncertainty  is  a  com¬ 
bination  of  uncertainties  in  tunnel  test  conditions,  model  positioning,  and 
model  instrumentation.  Since  the  lift  coefficient  value  for  the  flight  con¬ 
dition  of  interest  rarely  appears  as  a  test  point,  a  fitted  curve  must  be 
constructed  from  the  data.  This  curve  fit  process  is  another  source  of  possible 
uncertainty.  This  report  investigates  the  uncertainty  associated  with  wind 
tunnel  drag  data  and  presents  a  method  that  can  be  used  to  calculate  it  for  a 
specific  flight  condition  of  interest. 
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SECTION  II 


GENERAL  METHOD  OF  UNCERTAINTY  CALCULATION 

Uncertainty  associated  with  CD  at  a  given  condition  is  caused  by 
uncertainties  in  the  measurement  process  and  those  induced  by  the  curve  fit 
procedure.  This  requires  that  the  uncertainties  associated  with  each  factor 
be  propagated  to  arrive  at  a  final  result.  The  method  chosen  here  is  the  Taylor 
series  method  of  error  propagation. 

A  derivation  of  the  Taylor  series  method  can  be  found  in  Reference  1.  Some 
of  the  assumptions  used  in  the  derivation  are: 

1.  Response,  Z,  is  defined  as  a  function  of  the  measured  variables 

xi»  x2,  ...»  xn. 

2.  Z  is  continuous  in  the  neighborhood  of  u  .  y„  .....  y„  .  y„  , 

Xj  X2  Xn  Xj 

u  ,  . . . ,  yv  are  the  mean  values  associated  with  Xj,  x2,  ...,  x  which  all  have 
x2  xn  n 

error  distributions  about  the  point  of  interest. 

3.  Z  has  continuous  partial  derivatives  in  the  vicinity  of  u  ,  y  ,  .... 


4.  xj,  X2»  ....  x  are  independent  of  each  other. 

5.  (uv  -  xj,  (y  -  x2),  ....  (u  -  x  )  are  small  or 

xi  x2  xn  n 

Vz  ,  a2z  ,  ....  32z  are  small  or  zero. 

- ?  - 5“  - 2 

3Xj  3X2  3Xn 

The  assumptions  will  be  satisfied  if  the  functions  considered  are  restricted 
to  smooth  curves  near  the  point  of  interest  with  no  discontinuities  and  with  higher 
order  derivatives  either  small  or  zero. 

The  results  of  the  derivation  show  that: 


V 


If 

z  =  f  (xlt  x2,  ....  xn)  (1) 

then 


S(Z)  ='  {  HMb  S(Xl)  ]  "+[ 


3  Z 


3  X 


SX2  sl»2>)2  +  •••  + 1 H;  S(*n)  1 2 ) 


(2) 


where  S(Z),  S( X2 ) ,  S(X2),  ...,  S(X  )  are  the  precision  indices  of  the  response, 
Z,  and  the  variables  Xj,  X2,  Xn>  The  precision  index  is  the  computed 

standard  deviation  of  the  measurements  (i.e.,  random  error).  It  is  defined  as 

1  "  '  (3) 


S  = 


iw  ^  (,1i  -  *>2 


X^  and  x  are,  of  course,  the  value  of  x  of  a  particular  point  and  the  mean 
value  of  x,  respectively.  If  the  sample  size  is  large,  the  precision  index  is 
approximately  equal  to  the  actual  population  standard  deviation  (a)  associated 
with  the  random  variable  Z. 

Depending  on  the  confidence  level  attached  to  the  response,  the  uncer¬ 
tainty  is  simply  a  function  of  the  precision  index  or  standard  deviation. 

This  assumes  that  no  bias  or  systematic  errors  are  present.  For  drag  increments 
caused  by  external  stores,  this  is  a  reasonable  assumption.  The  drag  increments 
are  calculated  as  differences  between  data  taken  during  one  test  with  the  same 
model  and  instrumentation.  While  bias  errors  may  be  present,  they  are  approximately 

equal  for  different  model  configurations.  When  the  increment  is  determined,  the 
bias  errors  should  drop  out.  This  assumption  is  of  critical  importance  in  the 
following  discussions.  It  would  not  be  true  if  raw  data  from  separate  tests 
were  compared  or  if  different  instrumentation  was  used  during  one  test. 

For  the  case  of  interest,  Cl  and  Cq  are  both  random  variables  since  both 
are  measured  during  testing.  With  the  assumption  that  measured  values  of  Cl 


associated  with  C[j  on  the  fitted  Cl  -  C^  curve  is  found  as  follows. 
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Figure  1  shows  a  number  of  data  points.  From  aerodynamic  considerations. 


Cq  is  considered  a  second  order  polynomial  function  of  C|_.  That  is: 

CD  =  ao  +  3j  CL  +  32  CL2 

Using  the  regression  coefficients  aQ,  aj,  and  a2,  the  values  of  Cp  at  the  data 
point  ( C(_  1 )  and  the  flight  condition  of  interest  (C|_  )  can  be  found. 

%  '  ao  *  a.  CL,  +  *2  151 

cd2  ■  a0  +  a>  cl2  *  a*  cl2  (6) 

The  difference  between  the  values  of  Cp  is: 


CD2  "  %  ACD  3l  <CL2  "  CL^  +  32^CL2  "  CLi  '  ^ 

C[)?  can  be  redefined  in  terms  of  Cqj  and  ACp. 

C02  ■  CD,t4CD  =  C0,  +  a>  (CL2  -  CL,>  *  a2  <Cl2  -  Cl(>  <8> 

Now  assume  that 

V  f  <V  V  V  <9) 

Using  the  Taylor  series  method  of  error  propagation: 

?CD  ?Cq2  3CD2  2  J-2  (10) 

s(cn  )  -•  {[ - -  S(C  )  ]  +5 - S(C,  )  ]  +[—  S(C ,)]  } 

D2  uac  °l  1  3  C.  Ll  3  C.  12 

Dj  Lj  L 2 

The  partial  derivatives  in  Equation  (10)  can  easily  be  evaluated  from  the 
definition  of  Cn  in  Equation  (8). 


They  are: 


9CD2 

"3C7 


=  1 


K. 


9c, 


-  (aj  +  2a2  Cp  ) 


(11a) 


(lib) 
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Z  =  «!  +  2a2  C. 

l2 


(lie) 


The  terms  S(Cn  ),  S(C,  ),  and  S(C,  )  in  Equation  (10)  require  interpretation. 

Uj  Lj  l2 

They  are  the  precision  indices  of  Cq  and  CL  at  point  1  and  that  of  CL  at  point 
2.  The  at  point  2  is  simply  the  desired  flight  condition.  It  can  be  thought 
of  as  a  mathematical  and  not  a  random  variable  and  can  be  exactly  defined. 
Therefore: 

S(C,  )  =  0  (12) 

*-2 

C.  is  a  measured  data  point,  and  because  of  this,  it  is  a  random  variable 

Li 

with  a  mean  and  standard  deviation.  S(C^)  can  be  thought  of  as  the  uncertainty 
associated  with  the  measurement  process.  If  S(CL)  is  constant  or  it  changes 
slowly  in  the  vicinity  of  the  flight  condition  of  interest,  S(C.  )  should  be  a 

Li 

close  approximation  to  the  value  that  S(C. )  would  have  if  C.  were  a  measured 

l  l2 

data  point.  Figure  2  indicates  that  S(C^)  does  not  vary  much  over  the  angle  of 

attack  range  of  interest  (a  =  2  -  6°).  Because  of  this,  S(C.  ),  where  C,  is 

Li  Li 

the  data  point  nearest  the  flight  condition  of  interest,  can  be  used  to  indicate 
the  uncertainty  associated  with  the  measurement  process.  Calculation  of  S(C^) 
is  discussed  in  Section  III. 

Recall  that  the  value  of  Cn  used  in  Equation  (8)  to  define  Cn  was  not  the 

L>i  l>2 

measured  value  at  C,  ;  it  was  the  value  obtained  from  the  curve  fit  equation. 

Li 

In  other  words: 

Cn  =  an  +  a!  C.  +  a2  C. 2  (13) 

Dl  O  Lj  M 

2 

Therefore,  SfCp^  is  actually  the  precision  index  of  Cp ,  S  ,  introduced  by  the 
curve  fit  process.  A  method  for  calculating  S(C^)  is  found  in  Section  IV. 

Given  the  values  of  S(C.  )  and  S(Cn  ),  the  final  calculation  of  the  uncertainty 

Li  ui 

in  Cg  at  a  given  flight  condition  can  be  performed  as  follows: 
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u  <CD)  =  tn  S  (14) 

t  is  the  nth  percentile  point  for  the  two  tailed,  student  "t"  distribution. 
Various  percentiles  or  confidence  levels  can  be  chosen.  The  following  discussion 
is  based  on  a  95-percent  confidence  level.  The  ninety-fifth  percentile  point, 
to ^ 95 ,  depends  on  the  number  of  degrees  of  freedom,  which  is  a  measure  of  the 
sample  size  used  to  determine  the  precision  index.  For  S(C^),  t0ig5  is  usually 
taken  equal  to  2  since  S(C^)  is  based  on  a  large  sample  size.  This  value  is  not 
appropriate  to  use  for  S(Cp).  The  Aircraft  Compatibility  Branch  uses  a  five-point 
curve  fit  around  the  appropriate  C^  to  find  C^.  For  this  case,  the  degrees  of 
freedom  are: 

v  =  n  -  k  -  1  =  2  (15) 

where: 


n  -  number  of  points  used  =  5 

k  -  number  of  variables  in  regression  equation  =  2  (i.e.,  C^,  C^2) 

For  two  degrees  of  freedom: 

^0.95  =  4.303 

The  equation  for  uncertainty  can  be  rewritten  as: 

U(CD)  -  Uto.95  S(CD  )]2  +  [to. 95  (a}  +  2a2  CL  )  (16a) 

CD  1  CL  1 

=  {4.3032  S(Cn  )2  +  22[(a  +  2a  C.  )  S(C.  )]2}h  (16b) 

uj  1  2  L 1  *-l 

The  value  U( C^)  is  the  uncertainty  in  drag  coefficient  at  a  given  flight 
condition  for  a  given  configuration.  More  important  is  the  uncertainty  in  the  drag 
increment  between  the  clean  aircraft  and  the  aircraft  with  the  external  stores. 

The  increment  is  defined  as: 


ACr 


'Dr 


-  cr 


Stores  Clean 
From  Reference  1,  the  uncertainty  is: 

)2  +  U(Cr 


u(acd)  =  {U(CD 


Stores 


Clean 


)2> 


(17) 


08) 
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SECTION  III 


DATA  POINT  PRECISION  INDEX 

Section  II  described  a  method  to  calculate  drag  data  uncertainties  for 
a  point  on  a  fitted  C^-C^  curve.  Because  of  the  random  nature  of  the  assumed 
independent  variable,  C^,  the  final  calculation  for  U(Cg)  requires  a  value 
for  the  precision  index  of  C^.  This  can  be  done  in  the  following  manner: 
can  be  defined  as: 

-^1  [FN  (cos  a  cosIM  +  sin  a  sinIM)  (19) 

+  FA  (cos  a  sinIM  -  sin  a  cosIM)  ] 

If 

CL  =  f ( X] ,  x2 ,  . . . ,  xn)  (20) 

then  the  Taylor  series  method  of  error  propagation  yields: 

S(C.  )  =  {[  S(Xl)p  +  [^L  S(x2 )  ]2  +  ...  +[  S(x  l]2)5*  (21) 

L  8Xi  3x2  3Xn 

It  is  apparent  that  the  equations  become  more  and  more  complex  as  the 
number  of  independent  parameters  increases.  In  addition,  these  variables 
must  be  independent  or  nearly  so  of  each  other  for  Equation  (21)  to  hold. 
Assumptions  need  to  be  made  on  the  nature  of  the  measurements.  These  are: 

(1)  The  uncertainties  associated  with  some  of  the  independent 
parameters  are  sufficiently  small  that  their  effects  on  the  uncertainty 

of  the  dependent  parameter  are  negligible.  This  is  assumed  to  be  the  case 
for  IM  and  A. 

(2)  The  measured  forces,  FN  and  FA,  are  independent  of  one  another. 

Using  these  two  assumptions,  it  can  be  said: 

CL  =  f  (FN,  a,  FA,  Q)  (22) 

Therefore: 

S(C.)  =  {f^L  S(FN)  ]2  +  [^L  S(a)]2  +  [^L  S(FA)]2  +  [\  S(Q)]2}^ 

3FN  ?a  dFA  3  Q 
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After  evaluating  the  partial  derivatives.  Equation  (23)  may  be  rewritten 


as: 


S(CL)  =  ([(cos  a  cosIM  +  sin  a  sinIM)  ^jp-]2  +  [CQ  S(a)]2 

+  [(cos  a  sinIM  -  sin  a  cosIM)  ]2+  fCL  S(Q)j2}>s  (24) 

For  a  given  test  point,  a,  IM,  Q,  A,  Cp,  and  are  defined.  S(FN),  S(FA), 


S(a),  and  S(Q)  are  functions  of  test  conditions  and  instrumentation.  Equations 
for  these  parameters  can  be  developed  in  a  manner  similar  to  that  used  for 
Equation  (24).  S(a)  and  S(Q)  depend  on  the  wind  tunnel  where  the  data  are  taken. 
S(FN)  and  S(FA)  are  dependent  on  the  tunnel  and  model  instrumentation.  Equations 
to  calculate  the  quantities  for  a  specific  case  of  interest  should  be  obtained 
from  the  applicable  test  facility.  The  computer  program  described  in  Section 
V  uses  equations  that  apply  to  the  Aerodynanmic  Wind  Tunnel  (PWT/4T)  at  Arnold 
Engineering  Development  Center  (AEDC). 
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SECTION  IV 


CURVE  FIT  PRECISION  INDEX 

A  method  to  determine  the  precision  index  of  CL  for  use  in  Equation  (16) 
was  just  described  in  Section  III.  The  other  input  required  to  solve  for  U(Cp)  is 
SfCp^).  Recall  that  it  was  defined  as  the  precision  index  of  the  curve  fit  process. 
To  calculate  this  quantity,  certain  assumptions  on  the  nature  of  the  CQ  data 
are  made.  These  are: 

(1)  Cp  is  a  random  variable. 

(2)  At  a  given  C^,  possible  values  of  Cp  are  approximately  normally 
distributed. 

(3)  Over  the  CL  range  where  the  curve  fit  is  applied,  S(Cp)  is  nearly 
constant.  This  allows  a  simplified  curve  fit  procedure  to  be  used  (see  pp  106-108, 
Reference  2). 

Since  it  is  possible  to  calculate  S(Cp)  for  a  given  data  point  in  a  manner 
similar  to  that  used  for  S(CL),  the  validity  of  the  third  assumption  can  be 
verified.  Figure  3  shows  a  typical  example.  The  angle  of  attack  region  of 
interest  is  about  2  to  6  degrees.  As  the  figure  indicates,  S( Cp)  does  not  vary 
significantly  over  this  range. 

Recall  from  Section  II,  the  C^-Cp  curve  is  assumed  to  be  of  the  form 

CD  =a0  +  ai  CL  +  a2  CL2  (25) 

For  a  second  order  curve  fit,  the  normal  equations  to  solve  for  aQ,  ai, 
a2  are,  in  matrix  form: 

[B]  [a]  =  [g]  (26) 

where 
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1.1  Di 


[g]  = 


1*1  %  C[)i 


\  CL •  CD • 

1=1  1  l 


.  ''L. 

i=i  Li 


n  2 

•Z  CL . 
i=i  i 


*  CL. 

i=l  l 


n  2 

I  CL . 
1=1  l 


*  CL. 
1=1  l 


I  C. 
i=i  Li 


I*1  i 


E  c 

1=1  Li 


Solution  for  the  regression  coefficients  is  given  by: 

[a]  =  [B]  _1  [g  ] 

Where  [B]"1  is  the  inverse  of  [B] .  For  a  five-point  curve  fit,  n  is  equal  to  five. 

In  this  study,  the  value  of  S(Cq)  is  based  on  the  mean  response  ^  2* 

D  L ,  L 

The  term  ,r  r  may  be  thought  of  as  the  true  value  of  Cn  for  a  given  C.  . 
CD|CL,CL2 

In  other  words,  if  many  values  of  Cq  were  measured  at  a  given  C^,  the  mean  of 

C_  would  tend  to  yr  ,r  r  as  the  sample  size  increased  to  infinity.  Based 
D  cdICL,CL2 

on  this,  S(Cp)  is  defined  as: 

S(CD)  =  s  /[xo']  [Bp*  [xo  ]  (31) 
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where 


s*  estimated  sample  standard  deviation 
n 


E  fcD-  "  (ao  +  a*  ^Li  +  a2  cLi 
i=i  1 


2)]2 


n-k-1 


v-  n 


-FC 


■FC 


2] 


(32) 


(33) 


(34) 


lFC 
CL  2 

lfc 

The  FC  subscript  indicates  the  flight  condition  of  interest  where  the  calculation 
is  made.  For  a  more  detailed  discussion  of  the  method  used  to  calculate 
S(Cp),  see  a  statistics  text  such  as  Reference  3. 
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SECTION  V 

COMPUTER  PROGRAM  AND  SAMPLE  CALCULATION 

Appendix  A  contains  a  computer  program  developed  in  the  Aircraft  Compatibility 
Branch  to  calculate  uncertainties  in  drag  coefficient  where  Cp  is  found  from  a 
fitted  C^-CD  curve.  Required  inputs  are  explained  in  the  program.  The  uncertainty 
can  be  determined  for  confidence  levels  of  99,  95,  90,  and  80  percent,  i.e., 
x  percent  of  measured  data  points  should  lie  within  the  uncertainty  band  above 
and  below  the  fitted  curve.  The  confidence  level  is  changed  by  usinq  different 
t  values  as  explained  in  Section  II. 

As  noted  earlier,  the  equations  used  to  define  S(C^)  are  applicable  to  the 
four-foot  transonic  wind  tunnel  at  AEDC.  Changes  would  be  required  for  other 
tunnels. 

Assumptions  are  made  to  simplify  the  equations  for  S(C^).  These  include: 

(1)  Only  aircraft  pitch  excursions  are  considered,  i.e.,  3=0. 

(2)  Uncertainty  in  angle  of  attack  is  constant  and  equal  to  0.1  degree. 

(3)  S(a)  equals  one  half  the  uncertainty  in  a  (0.05  degree). 

(4)  A  small  sting  roll  angle  is  assumed  (0.04  degree). 

(5)  Model  weight  does  not  vary  between  configurations. 

(6)  To  a  first  approximation,  model  angle  of  attack  is  equal  to 
sting  pitch  angle. 

(7)  Balance  uncertainties  (S(FNB),  S( FAB ) ]  are  functions  of  normal  and 
axial  forces  only  (i.e.,  no  side  load  or  rolling  moment  interactions  are  present). 

(8)  There  is  no  model  roll  angle  relative  to  the  balance. 

(9)  The  precision  indices  of  model  weight  measured  by  axial  and  normal 
force  gages  are  equal  to  the  precision  indices  of  the  balance. 

(10)  Tunnel  total  pressure  is  less  than  1500  lbs/ft2. 

(11)  S(M)  is  constant  and  equal  to  0.002  which  is  40  percent  of  the 
uncertainty. 
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The  assuirptions  are  reasonable  for  the  data  of  interest.  Only  pitch  excursions 
are  considered  because  the  performance  problem  concerns  mainly  flight  at  a  constant 


angle  of  attack. 

Appendix  A  contains  sample  calculations  of  U(Cg)  for  a  specific  case.  The 
example  is  for  a  clean  aircraft  at  M  =  0.8  and  CL  =  0.3.  If  the  equation  for 
U(CD)  is  examined  using  the  values  of  S( Cp)  and  S(C^)  from  the  example, 
it  shows  that  most  of  the  uncertainty  in  Cp  is  a  result  of  the  curve  fit.  Note 
all  the  confidence  levels  are  shown  in  the  example,  and  U(Cp)  decreases  as  the 
confidence  level  is  lowered.  To  be  more  precise,  the  values  S(FN),  S(FA),  S(a), 
and  S(Q)  used  to  find  S(C^)  should  be  redefined  for  each  confidence  level 
instead  of  using  the  95-percent  values.  Since  the  contribution  of  S(CL)  to  U(Cp) 
is  not  large,  this  is  not  significant. 


SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

A  method  has  been  presented  that  shows  how  drag  coefficient  uncertainties  can 
be  calculated  for  wind  tunnel  data.  Because  of  the  need  to  examine  Cp  as 
a  function  of  Cp,  the  uncertainties  in  Cp  are  due  both  to  the  uncertainties  in 
Cp  and  those  of  the  curve  fit.  While  the  method  of  calculating  the  precision 
index  of  the  curve  fit  is  general,  that  used  to  determine  S(C^)  will  depend  on 
the  wind  tunnel  and  particular  test  instrumentation.  Consultation  with  the 
applicable  test  facility  will  be  necessary  to  work  out  suitable  equations  for 

s(cL). 

The  example  given  in  Appendix  A  indicates  that  the  uncertainty  can  be  fairly 
large  for  a  high  confidence  level.  Reducing  the  confidence  level  results  in  a 
considerably  smaller  uncertainty.  It  would  be  very  useful  if  data  uncertainties 
associated  with  flight  tests  were  given  a  thorough  analysis.  Determination  of 
these  uncertainties  would  increase  the  confidence  in  performance  estimates  based 
on  flight  test  data. 
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APPENDIX  A 


COMPUTER  PROGRAM  FOR  UNCERTAINTY  CALCULATION 
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it  W  -  •»;  v 


74/74 


CPT  =  C  T  PACE 


FTN  4.3+510 


1 


5 


n 


PROGRAM  UMCF  9(  ISO»jT,0\JTPUT,TAPFF=INP'jT,TtPr  6=CUT»UT) 


C  TH;s  P^GRAH  CALCULATE^  THE  UNCEcTAINTr  IN  fPAG  COEFFICIENT 
C  FOF  A  SPECIFIC  FLIGHT  CONDITION  OF  INTEREST  ON  A  FITTEO 

C  CL-CP  CU9Vr.  THE  UNCERTAINTY  IN  CO  IS  HA OE  JP  Of  A 

*  roHRTNATION  OF  UNCERTAINTY  IN  lift  COEFFICIENT  and  the 
C  IJNf  F  9  T  AI  N  T  Y  INTRODUCED  bf  CU^VF  FITTING  CL-CO  OATA. 

c 

C  the  following  constants  are  sro»n  in  diockoata. 

n 


lr' 


?F 


3  0 


C  PHII-ASSNHFO  STING  dOU  ANGLE  DUPING  A  PITCH  SWEEP  (DEG) 

C  spmi i-standapd  deviation  OF  STING  ~OLL  ANGLE  (DEGI 

C  SI-CONSTANT  DEPENDING  ON  PARTICULAR  WIND  TUNNEL 

C  w-AIPCPAFT  W  OflF  L  WEIGHT  ( L  9S  ) 

C  A^.Ao-CONSTANTS  USED  TO  OBTAIN  STANOARO  DEVIATION  OF  TUNNEL 

C  TOTAL  QPfSSUFE.  OEPENOS  ON  TUNNEL  OF  INTEREST, 

r  Sm-STANDARO  HE  V I A  T  I  f'N  OF  mACh  IMUMRFP 

C  SAtPHA-STANOARO  DEVIATION  CF  ANGLE  CF  ATTACK  (RAO) 

c  AtP-ANSLC  °  F  T  Wf  E  N  A  I  °OP  AF  T  WATER  LINE  AND  HALANCt  AXISIOEG) 

r  J-AIACPAFT  MODEL  REFERENCE  AREA  (SC/FT) 

C  A1 7, A16-CONSTANTS  USfO  TO  nE  T  ERH  lNr  STANDARO  DEVIATION  OF 

C  TUNNEL  STING  PITCH  ANGLE.  OEPENOS  ON  WINO  TUNNEL. 

C  Kl,XC,)(3,Y9-PALANCr  CONSTANTS 

C  A*1 .  A<15  « A K2R- BALANCE  CONSTANTS  ( C  A  LL  FD  Kl,<15.<29> 

C  SF AN2.CC AA2-OAL ANCF  AXIAL  Fo^rr  GAGE  UNCERTAINTIES  WITH 

C  APPLIED  NORMAL  AND  AXIAL  LOAD  Rt  SPEC  T I VEL  Y  • 

C  SFNN2tSFNA2-r,ALANCE  NCPHAL  FORCE  GAGE  UNCERTAINTIES  WITH 

C  APPLIED  NORMAL  AMO  AXIAL  FORCE. 

c 

c 

C  T  HE  FOLLOWING  VALUES  APF  CALCULATro  *-©0*  ThE  INPUT 


35 


O 


4  5 


r  0 

55 


r 

c 

c 

c 

c 

c 

u 

r 

c 

c 

c 

c 

c 


c 

c 

c 

C 

n 


cf LPI-STANOAPD  OEVIATION  CF  STING  PITCH  ANGLE 


JFA- 

SFK'- 

SQ- 

SCL  T  S- 


~  AXIAL  FORCE 

-  NORMAL  FORCE 

•  DYNAMIC  PRESSURE 

-  LIFT  COEFFICIENT 


T«e  FOLLOWING  INPUTS  APR  fqc?  THE  TEST  POTNT  CLOSEST  TO  THE 
flight  condition  D f  interest. 

ALPHA-AIPCF Aft  A  NCl  r  OF  ATTACK 
A-ACm-haCH  NJMRFP 

ot-tunnel  total  pressure 

0- T U NNFL  DYNAMIC  PRESSURE 
CLTS-LTFT  COEFFICIENT 
C"T$-ORAG  C OEFFICIEM 

°0  I-  V  AL  UE  of  LIFT  COrFFlCIENT  AT  THF  FLIGHT  CONDITION 

OF  interest. 


CONFI  IS  THf  c  CNF I  Or  NCE  LEVEL  SELECTED. 

VALUES  DF  .99,  .05,  .9,  AND  .*  ARE  AILOWEO. 

THf  FIVE  CL  AND  CD  VAL’JFS  ARE  THE  ONFS  USED  FQR  THE  CURVE  FIT 


GO 


F  5 


7  n 


75 


r.  XC.XOP.R.G,  AND  AC  APE  MATRICES  USED  TP  097  AIN  RCGRESSION 
0  COEFFICIENTS  AND  CALCULATE  THF  STANDARD  DEVIATION  DUE  TO  FIT 

C 

C  AQ,A1,A?  ARF  TH£  ©EGRESSION  COEFFICIENTS  CALCULATEO  FO*  THE 
C  F I  T  T r  P  CURVE  WMF  R£  IT  IS  OF  THE  FQRM 

C  CP  =  A 1  ♦  Al *CL  ♦  A2MCL"2) 

r 

C  SF I T 2  IS  Thc  SDUARE  OF  THE  STANDARD  DEVIATION  OUE  TO  CURVE  FIT 
C 
C 
r 

C0MHON/A/°HIIfS°HlI.SL*W*A5»A6.SMtSAL  PHA • AIM.A.A17.A10 
C0MH0N/rf/X9. CK15. $F AN2.SF AA2 
CDHiON/C/Xl,X2*X3tAKl,AK29*SFNN2,SFNA2 
COMMON/n/SALPr.SFA.SFN.SO^SCLTS 
COMMON/ E/ ALPHA* A M AC H.PT, O.CLTSfCOTS.POI 
CPMH0N/F/0,L<5>  .CD(5)  ,XC(3)  .X0P<3>  ,0(3,5>.G(3>,ACt3> 
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V 


Mom*  oncer 


TM/T it  0FT*0  TRACE 


FfN~*iO»tiR- 


COEFFICIENT /lH0»5F12.<*t 


COMM  OH/  5/  AO  ,  Ai  ,  A27SFITF - - 

COMM ON/ H/ AO ( 1I«AC(3I»AN(3>»Z(3I»X11J 

10  READ  <*,!>  ALRHA,ANACM«-PT»Q»CCTS«COTS»POI - 

1  FORMAM7FH.5) 

IF (EOF (51 1 99* 3 

3  REAO (5*81  COMF1  —  -  - 

8  FORMAT (F10*5I 

-  rail:?!  ISk!«*hl) - 

2  FORMAT  (5F10 • 51 

XA1L  SOCL  - 

CALL  SOFIT 

WRIT  E(6« 111  AMACH 

11  FORMAT  <  1*1,1 3MMACR  ROMBER  *oFir«vOT - 

WRITE C6* 12) 

u  5s?wu;-«r«i!?i"?.iw£i  MTAi«’,J-SL - 

13  FORMAT  C/SX • 26HVALUES  OF  LIFT  COEFFICI ENT/lM0t5F12# 4) 

WRITE  (6*141  _ - _ _ 

14  FORMAT (/5X.26HVALUES  OF  ORAG  COEFFICIEMT/1HO, 5F12 . 41 
WRITEI6* 15 1  POI 

15  FORMAT  IF //Sit rtttttff  T  OOEFFIOfFRF-OF  IRTERES T  ^»F10.41 
WRITE  C6»  71  A0« A&  « A2 

T3t^(4:>a3!l?l!»K6?? 

WRITE(6«41  SCLTS 

_ FORMAT (//IX. 24KSTAMPAR0  DEVIATION  IM  CL  g<F12«4l _ 

SFIT*SQRT(SFIT2) 

WRITE (6. 5 (  SFIT  _ 

—9"  FORMAT (//IXrgTHSTRHtWHr  DEVIATION  ~IN  FIT  »,Ft2.M - 

conf*6onfi*ioo. 

If  1^1(llt//!i*i1?MC5?iFTBEMCC  tMTERVAL  ~*,P5.S,2H  TJ - 


ID  ru^TRi  9  vCiriuunr  iucnuc.  iniCKVHi  *tr»«B»cn 

. . IHESSi-ffcfcUJJS  &il! _ 

IF(CONFI.EQ.0.9>  CO  TO  182 
IF (COMF1 >£0* 0,8)  SO  TO  1D3 

*•* tes?:?!’ - - 

- 

TWAL 2*2* 

_ cn  fa  9aa  _ _ _ _ _ _ _ _ 

102  TVAL1*2.92 

TVAl 2*1#65 _ 

113  ?TalT*iT1#6 

"Till — varV»(Hal1  - 

z  VAR2«(t}al|**1.I*((A1t2.*A2*CLTSI**2.»*(SCLTS**2.» 

- UCO»$ORT4WR4»VAR2I  - - 

MRITE(6,6)  UC0 

6  FORMAT (///IX, 3TMTHE  UNCERTAINTY  IN  ORAG  COEFFICIENT 

9ri?oTrt* 

_ £NIl- _ 


.F12.5. 


«,F12»RI 
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SUBROUTINE  SOCl 

-THIS  SUBROUTINE  CALCULATES-SCCTS.  -THE  STRNOARO  DEVIATION 
LIFT  COEFFICIENT  FOR  THE  DATA  POINT  NEAREST  THE  FLIGHT 


CONDITION  OF  I nterest. 


_ 4sag8a^^ss?i&f 6  -  s,^5^r^  ” 

COMHON/C/X1, Xg#X3, AK1*AK29*SFNN2,SFNA2 
COMHON/O/SALPI «SFA«SFNiSCUSCL  TS 

- eOMMW€/AlPHAtAMACH%PTvfrret«4CQfS«P0i  - 

SALPI*A17*A10*ALPMA 

SFAB=.5*< (SFAN2**2.^SFAA2**2.)**.5) 

- swisraar  — - — - - - — 


A17, A1S 


SFAG2=(SL*X9*AK15I**2.»SFAB**2. 

-  . -SF ASJ2=  W*COSO  (ALPHA) -«$AL*XM  t£  ^USJMO  C ALPHAI*SIIAA*A2* 
SFA=  CSFAG2*SFAST2>***5 
SALPHA=. 000673 

- Sf^e^^^tfSFWNz^p.^sFH^e^e.^s^  - 

SFNG 2*<<SL/X1I**2.I M (X2*AK1) **2 . ♦ f X3*AK29>  *»2*> 

1  ♦SFNB**2. 

5FNB - - - - - - 

VAR1 N*  (W*SlNO(Al  PH A ) ^ CO SO  CP HI I ) •SALPI )**2# 

- VAR2H=  ICOSC  (ALPHA)  *CQSDXPiiIIl*  _ _ 

VAR3N*CW*COSOf  AlPMA)*SINO(PHIII*SPHin**2. 

SFMS72- VAR1N*  VAR2N^  VAR3N 

- SfN*<SFNfrg+SPttSTgV*»^g — - - 

VAR1PS  (  f  I2**0/AHACH)M1.-1*2*  (ANACH»*2,)))/ 

1  <1-*.2MAMACH**2.  )))*SM 

— ^w-nrS5*w*FT - - - 

VAR2P* (Q/PTI *SPT 

Sq*  IVAR1P-* »2 . »  VAR2P •  - 

VAR1T® ( (COSO (ALPHA) *COSD)AIM) ♦StNO (ALPHA) #SIN0( AIMI I 
1  *SFN/(Q*A))**2, 

--  YAP?  T*  <$OTS*SAt PHAf  - - 

VAR3T* ( I  COSO C ALPHA) *STNO( AIM) -SINO ( ALPHA) *COSO ( A IMI) 


SCLtS*5 
RE TORN 
EMO 


S*SQRT(tfRRl7 ♦*AR2mAR37*VA**T) 


SUBROUTINE  SOFIT 


7 4/74  OPTs®  TRACE 


FT*  A.A*5iR— 


10 


15 


20 


2* 


SR 


R5— 


AO 


A5 


50 


55 


C 

C 

c 

c 

c 


SUBROUTINE  SDFIT 

THIS  SUBROUTINE  CAttBLATFS  T*F  RFORESSIU*  "COEFFICIENTS  ~TO~FTT~ 
THE  CL-CO  OATA  ANO  DETERMINES  THE  STANDARD  DEVIATION  OF 
.THE  CURVE  . FITj»  .  _  _ _ 

COHMON/E/ ALPHA*AMACH*PT»QjCLTS«C  DTS«POI 

COMMON/F/Cl ( 51 ♦CD  <  5  I ♦ XO  C3J  ♦XRP4R1**44URUC434*AC4R1 - 


__C _ 


COMNON/G/AO*  Al« A2«SFIT2 

COMMON/H/AOt II tAL<3>* AMC3I  •ZI3I«X(1) 

- 


SUM21 


10 


20 


c 

c 

i 


SUmS1=0» 

sum*  i=o. 

OO  10  1*1*5 
SUM11=SUN11»CL (1 » 
SUN21=SUN21»CL  <11 **2. 
SUN31=SUN31*CLtII**3. 

SUM'*1  =  SUM<*1»CL  (!)•*!*. 

CONTINUE  _ 

Btt* ll 55 • 

BU,2I=SUMU 
BU*3I*SUM21 
BI2, II *SUM11 
812, 2)=SUN21 
Bt2,BI=$UH3t 
(M3,1I*SUM21 

915»l!=S^" 

B<3.3I  =  SUI< 

SUM12*0, 

SUH22=0. 

SUNS  2=0, 

00  20  1*1,5 
StfN12*SUNi**CD<II 
SUM22=SUM22*CL II I  *CO (II 

G(  II  =SUN12 

G  (21 -SUM22  . -  - 

GI3I  =SUM32 

CALL  MINV(B«3,AD,AL*AMI 


CALL  TO  SYSTEM  ROUTINE  TO  INVERT  MATRIX  B 


CALL  GmPRB  (6~,G, ATT,3,3,‘i) 


MULTIPLY -INVERSE  OF  9  1Y  6-  EG- OBTAIN  MATRIX  -AC- 
A0  =  ACUI 
A1*AC(2I 

A2=AYJtS^  - 

SUM3=0. 


60 


SunlisfaK^O  m  -A 0 -Ai*CL  IIJ-A2*  CCL  1 

1  **2. 

CONTINUE  - - 

SSD*SUM3/2« 


XOPt 11*1# 

TRPtF1*P or 

XOI  - 


mw 


71 


XflPC  31 *P3I 

M.W  ‘ 


AO  CONTINU 


C 

* 


RO I X  OR  ♦ 


£a° 

MULTIPLY  XQP  BY  INVERSE  OF  B  TO  OBTAIN  Z 
CALL  GMPRO(Z,XO,X,1,3,1I 


— 7$ 


5CTIPLY”2RY  XO~T0B STAIN  X.  X  IS  A  ONE  COLUMN  ONEROU 
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60 


suBROOTine-  sgfit— 


7  «r77*  f**GE 


TT - DFVTSTTO>r  DUE  TO  CURVE  TIT." 

C 

— Xis*  III  - 

SPIT  2=SS0*Xi 
RETURN 


24 


St  OCK  OATABtKOAf. 


74/74  — 


WH  »«-«♦» i»- 


“BlOCKOm . . 

COHMON/C/X1 , X2»  X  3. AK1»AK29,SF  NN2.SFNA2 
DATA  PHI I, SPHII.SL  »H. A5.A6.SN.SALPHA. AIH, A,  I 


10 


t  '  •  0 4,  .04, 1  •  6 •  33" i«  •!'  . 
1gATASX^AKlS|SF^M2^P^A^^|.^3g|g^«^|/ 

1  .393,^39^4.1 89. 176/ 

-&W>  - - 
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V 


WIND  TUNNEL  DATA  POINTS 
”  VALUE S~OF”  LIFT  COEFFICIENT 

- .3231 —  »1?T« - 

_ VALUES  OF  DRAG  COEFFICIENT 

.0152  .0169 

--  LIFT-ODEFF10IENT  OF  ENTE*EST  « 


CURVE  FIT  COEFFICIENTS 

AO”  *  .  .01684 

A1  *  -.06606 

A2  *  — -  .225  97 

ST  ARUARO  DEVIATION  rNHCL  =  '* 

STANDARD  DEVIATION  IN  FIT  =  ~ 

CONFIDENCE  INTERVAL  =  <»9.  X 


3324 - .4925 - .6374 

0272  .0466  .0603 


*3905 - 


TOT!  37" 

.0014 


THE  UNCERTAINTY  IN  OR  AH  COEFFICIENT 


.9139 _ 


ffnvn  Nunut*  * 


.799 


NINO  TUNNEL  DATA  POINTS 

VALUES"  OF  LIFT  COE F PICIFN T 

- .0*31 -  .IFF*  '  *332*  - *4*25 

_ VALUES  OF  DRAG  COEFFICIENT _ 

•0152  .0169  .0272  .0466 


- LIFT  COEFFICIENT  OF  I+JTEPEST  = - .3-003 


CURVE  FIT  COEFFICIENTS 

AO  _=  .01664 

A1  *  -.04606 

- a*  - ftmr —  —  - 

3TANUXRTT C EVI ATTW  T*T  CT~=*  YDOTJ 

STANDARD  0 EVIA T ION  IN  F I T  =  .0014 

CONFIDENCE  INTERVAL  =  95.  7. 


_ .4044 


-  THE  UNCERTAINTY  IN  -DRAG  COEFFICIENT  = 


NA04  NUH8£R-g - rTO* 


HIND  TUNNEL  DATA  POINTS 

values  'of* Lift'  coefficient 

- •  8  231 - .*770 - t 3324 - .4925 - iW» 

VALUE$_OF  DRAG  COEFFICIENT _ 

•0152  .0165  .0272  .0460  .0803 


LIFT  COEFFICIENT  OF— 2NTEREST-  * 


.3000- 


CUPVE  FIT  COEFFICIENTS 


AO  = 
A1  = 
AO— * 


.01684 
-• 04606 
“•-2200T- 


jnrwonnj  deviation  tn  ci  = 


STANDARD  DEVIATION  IN  FIT  = 


T0P73 

.0014 


CONFIOENCE  INTERVAL  *  90 .  X 


THE  UNCERTAINTY-^.  OR1&  -COEFF ICI EH-T s. -  .0041 


28 


HfrCN  NUMBER 


rm 


KINO  TUNNEL  OATA  POINTS 

VALUES  "oFliFT  COEFF IC I  ENT 

-  — <,231 - rtTTd - .3324 

_  VALUES  OF  DRAG  COEFFICIENT _ 

•  0 152  .0169  .0272 


Iff  OOfFFlOIf KT  OF  I  NTFREST  ~ -  .3000 - 


CURVE 

FIT 

COEFFICIENTS 

AO 

= 

•01664 

ii 

- s - 

-.04636 
-  .22597 

- STANDARD'  DEVIATION  TTTCt  =  .DU 3TJ 


STANDARD  0EVIATI0N  IN  FIT  =  .0014 


CONFIDENCE  INTERVAL  =  60.  X 


- TNE  UM&ERTA-IMT-T  IN  DRAG  COEFFICIENT  - - .04127 


.0466 


29 

(The  reverse  of  this  page  is  blank) 


